(N 



O 
(N 

H 



O 






> 



(N 



% 



Experimental verification of Boltzmann equilibrium 
for negative ions in weakly collisional electronegative 
plasmas 

Young-chul Ghim^ and Noah Hershkowitz^ 



^Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, 
g : 0X1 3NP, United Kingdom 

^Department of Engineering Physics, University of Wisconsin-Madison, 
CN ! Madison, Wisconsin, 53706, USA 

Abstract. Weakly collisional Ar-02 electronegative plasmas are investigated in a dc multi- 

dipole chamber. An electronegative core and an electropositive halo are observed. The density 
j^ . ratio of negative ions to electrons (a) in the nondrifting bulk is found to be 0.43. The profile of 

p ' a is found using both the phase velocity of ion acoustic waves and the drift velocity of positive 

j/5 ! ions determined by laser induced fluorescence. The experiment shows that negative ions are 

in Boltzmann equilibrium with a temperature of 0.06 ± 0.02 eV. Double layers are not found 

separating the electronegative core and the electropositive halo. 



Experiments have demonstrated that electronegative plasmas are separated into three re- 
gions yj-y]. These three regions are an electronegative core, an electropositive halo, and an ion 
lO ■ rich sheath region where the core and the halo in some parameter regimes Jol-lsll are predicted 

to be separated by a double layer. Braithwaite and Allen [9J, and Franklin [llOn argued that the 
lO ■ appropriate parameters are T^/T^ > 10 and ao = n_o/neo > (Te/T_)^/^. Here, Tg, T_, rieo, 

cn ! and n_o denote the temperature of electron and negative ions in eV, and the number density 

^ I of electron and negative ions in the bulk plasma, respectively. Berezhnoj et al. yD determined 

negative ion density profiles using photodetachment and Langmuir probe data in oxygen capac- 
itively coupled rf plasmas. They remarked that negative ions were in Boltzmann equilibrium 
for low pressure in comparing numerically calculated fluxes of negative ions and positive ions. 
Plihon et al. ^ reported the existence of double layers in electronegative Ar-SFg plasmas. 
However, these double layers were associated with a cusp magnetic field. The double layers 
with potential steps of 8 V^ for negative ion temperatures assumed to be at room temperature 
(~ 0.025ey) were not observed when electropositive plasmas were studied in the same system. 
Several authors have argued that it is appropriate to describe negative ions in low pressure 



electronegative plasmas by Boltzmann equilibrium [|llNl3ll. However, direct measurements 
have not been carried out to verify if negative ions are in Boltzmann equilibrium in electroneg- 
ative plasmas. Spatial profiles of negative ion concentration and plasma potential are difficult to 
determine accurately. This letter is the first experimental report to verify that negative ions are 
in Boltzmann equilibrium and that the double layer does not exist in weakly collisional elec- 
tronegative plasmas for low negativity {a^ = 0.43). The negative ion temperature is measured 



in the electronegative plasmas. 

The dispersion relation of ion acoustic waves (lAW) in electronegative plasmas can be de- 
rived from the general electrostatic wave dispersion relation [Il4ll . 
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where u is the frequency, k is the wavenumber, v^ is the drift velocity, and Vt is the thermal 
velocity. The subscript j denotes the jth species. cOp is the plasma frequency, and the sum 



is over all species. In electronegative plasmas with non-drifting ions, D'Angelo et al. lllSn 
predicted that there exist two modes (fast mode and slow mode), and two modes were observed 
by Sato et al. u6\ . If positive ions drift in electronegative plasmas, then each mode has two 
different phase velocities corresponding to the lAW traveling in the direction (parallel mode) 
and in the opposite direction (antiparallel mode) of the ion drift. 

The fast mode lAW dispersion relation can be derived from Eq. ([T]) by assuming tu/k ^ 
Vt+, u/k :$> Vt-, uj/k <^ vte, and kX^e <^ 1 where Ad is the Debye length. Here, the 
subscripts +, — , and e denote positive ions, negative ions, and electrons, respectively. Making 
these assumptions, the fast mode lAW phase velocities (Vphp) with nonzero Vd+ are 
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where cs = ^/Te/m^, and a = n^/nf.. 

The slow mode lAW dispersion relation can also be derived from Eq. ([T]) by assuming 
uj/k ^ vt+, ijjjk <^ Vt-, uj/k <C Vte kXoe ^ 1, and kXo- ^ 1- With the above assumptions 
and nonzero Vd+, the slow mode lAW phase velocities (Vphs) are 



VphS = Vd+±Cs, 
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In RHS of Eq. dS]), the '+' sign corresponds to the parallel mode, and the '— ' sign corresponds 
to the antiparallel mode of the slow mode lAW. 

The lAW dispersion relation in electropositive plasmas is derived from Eq. (d) by assuming 

u/k ^ Vt+, u/k <^ Vte, and kXjje ^ 1 with nonzero Vd+: 



Vph = Vd+ ± Cs, 



(4) 



where Vph is the phase velocity of the lAW in electropositive plasmas. The parallel (-I-) mode 
and the antiparallel (-) mode exist due to drifting ions. 

In electropositive plasmas, the location of the sheath edge is found to be where the parallel 
mode of the Vph is twice cs [Il7ll . It is equivalent to argue that the antiparallel mode of the Vph 
is zero at the sheath edge. The same argument has employed to locate the inner sheath edge 
of electronegative plasmas. Braithwaite and Allen \^ argued that at the inner sheath edge, the 
parallel mode of Vphs is twice C5a/(1 + «)/(! + aTe/T—), whereas the antiparallel mode of 
VphS becomes zero. This argument is valid if a double layer exists at the inner sheath edge. 

The experiment was carried in a cylindrical dc multidipole chamber 70 cm in length and 
60 cm in diameter (interior dimensions). A detailed description of the chamber can be found 
elsewhere [llSll . A stainless steel plate diameter of 15 cm. is biased to —30 V" in the center of the 
chamber to act as a boundary. Electronegative plasmas were generated by primary electrons 
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Figure 1. (Color online) Profiles of v^^ and (j)piasma from z = 15.0 cm to 0.2 cm. 



emitted from thoriated hot filaments biased at —60 V and the discharge current 0.8 A with em- 
ployed neutral gases of Ar at 0.10 niTorr and O2 at 0.30 mTorr. A planar Langmuir probe 
(LP) with radius of 0.3 cm was utilized to provide the electron number density and the effective 
electron temperature in the bulk plasmas. In this experiment, the LP trace showed the existence 
of hot and cold electrons. Thus, the effective electron temperature was used to find the Bohm 
speed, cs n\% . From the LP trace it is found that rig ~ 3.8 x 10^ cm~^ and Tg ~ 0.69 eV which 



results in cs ~ 1290 ± 50 m/ s with m+ being the mass of Ar. Use of the mass of Ar for m+ is 
justified in the following paragraph. 

A tunable diode laser was used to determine the profile of Ar+ drift velocity (fd+). The 
laser of 668.614 nm (in vacuum) excited the metastable Ar+ (3(i^F7/2) to Ap^D^i2 level, which 
in turn emits fluorescence radiation with a wave leng th of 442.72nm (in air). A detailed descrip- 
tion of the LIF system can be found elsewhere ll20ll . It is assumed that the metastable Ar+ is in 
thermal equilibrium with the ground state Ar"*" [12 ill . By scanning the wavelength of the laser, 
the ion velocity distribution functions (ivdf) were obtained from the Doppler shift. The values 
of Vd+ are found by taking the first moment of the ivdfs. Figure [T] shows the profiles of the Vd+ 
and the plasma potential {(ppiasma) measured by an emissive probe using the inflection point 
method in the limit of zero emission [12211 . The sheath edge was located at z = 0.6 cm. Here, z 
is the distance from the boundary to the diagnosed point. The measured Vd+ at z = 0.6 cm was 
1230 ± 70 m/s. This velocity agrees well with the calculated cs with the assumption of m+ 
equal to the mass of Ar. If there were non-negligible 0+ or 0^+, then the measured Vd+ would 
have been faster than the calculated cs 1231]. 

The values of the phase velocity of the lAW were measured using CW (continuous wave) 
excitation at 100 kHz with peak-to-peak voltage of 3.0 V^ and offset voltage of 0.5 V. A detailed 
description of the CW lAW setup can be found elsewhere [llVll . The measured phase velocity 
is denoted by VphM. The excitation was launched at the center of the chamber with the launch- 
ing grid diameter of 10 cm.. The LP is located between the launcher and the boundary which 
were separated by 20 cm.. The LP was positively biased so that it collected electron saturation 
current and was moved from z = 15.0 cm, to 0.5 cm,. Here, z is, again, the distance from the 
boundary. The signal collected by the LP was averaged by a boxcar averager with 5 kHz trig- 
gering, 30 nsec window width, and delay time varied from 0.0 to 10.0 /isec with the step size 
of 0.5 fisec. Figure |2] shows the profile of VphM from z = 15.0 cm, to 1.0 cm.. 

The electronegative core and electropositive halo regions are identified by comparing Vp^M 
to Eq. dU. In the electronegative core, the values of Vp^M are not equal to Vp^ in Eq. dU, whereas 
in the electropositive halo, these two velocities are equal. Figure |2] shows that these two ve- 
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Figure 2. (Color online) The profiles of Vp^M and cs + fd+ = Vph- The values of a are calculated 
using Eq. ^ with the measured VphM, Vd+ and cs- 



££ 0.1 

3 

E 

M 

n 

.^-0.1 



N,-0.2 
1-0.3 



----Linear Fit 



1 



^4. 









-3 -2 -1 

Log.(n(z)/n.o) 



Figure 3. (Color online) Negative ions follow the Boltzmann relation with T_ 



0.06 ± 0.02 eF. 



locities become the same at z = 2.0 cm. Thus, the electropositive halo starts at z = 2.0 cm. 
Consequently, from z = 15.0 to 2.0 cm., the measured phase velocity VphM corresponds to the 
fast mode of the lAW (Vphp) because this region is the electronegative core. Using Eq. dH) 
and having the measured values of Vphju (= Vphp), cs, and Vd^, a can be calculated. Here, the 
dominant negative ions are assumed to be 0~ (i.e., m_ = 16) [lli|3|]. Figure |2] shows the a 
profile and the existence of the electronegative core and the electropositive halo. 

To verify the Boltzmann equilibrium for the negative ions, the Boltzmann relation is exam- 
ined over the measured values of a and plasma potentials. Figure |3] shows that the negative ions 
agree with the Boltzmann relation with T__ = 0.06 ± 0.02 eV within the experimental errors. 
In addition, the plasma potential profile shows that there is approximately a 0.5 ±0.1 V drop 
from z = 15.0 to 2.0 cm,. This potential drop is sufficient enough to confine the negative ions 
with the temperature of 0.06 ± 0.02 eV within the core. 

An electronegative core and an electropositive halo were detected. Nevertheless, the plasma 
potential profile shows that there are no double layers. The fact that the negative ion fraction 
is gradually decreasing (Fig. |2]) rather than making an abrupt reduction suggests nonexistence 
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Figure 4. (Color online) RHS of Eq. ([S]) as a function of a with Tg 
0.08 eV (lines), and measured Vd+ as a function of a (square dots). 
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of double layers. The absence of double layers can also be verified by considering the positive 
ion drift velocity Vd+. Eq. ([5]) is the theoretical prediction of the Vd+ at the inner sheath edge 
when there is a double layer. 
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Comparing the measured v^^ with the RHS of Eq. ([5]) which is the Vp^s with nondrifting ions, 
Fig. m shows that Vd+ never reaches cs\/{l + «)/(! + aT^jTJ). Thus, it can be concluded 
that the double layer for the investigated condition does not exist. This is consistent with the 
theoretical prediction [|9lll0ll. 

In conclusions, using the measured drift velocity of positive ions and the phase velocity 
of the lAW in electronegative plasmas, the profile of the negative ion fraction is found in a 
dc multidipole chamber. The location where the electronegative core and electropositive halo 
are separated is found at z = 2.0 cm, while the electropositive halo terminates at z = 0.6 cm 
followed by the sheath. It is found that the negative ions follow the Boltzmann relation with 
a temperature of 0.06 ± 0.02 eV and that a double layer does not exist when negative ions are 
confined in the core of the plasmas. 
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